Evaporation of the pancake-vortex lattice in weakly-coupled layered superconductors 
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We calculate the melting line of the pancake-vortex system in a layered superconductor, interpo- 
lating between two-dimensional (2D) melting at high fields and the zero-field limit of single-stack 
evaporation. Long-range interactions between pancake vortices in different layers permit a mean- 
field approach, the "substrate model" , where each 2D crystal fluctuates in a substrate potential due 
to the vortices in other layers. We find the thermal stability limit of the 3D solid, and compare the 
free energy to a 2D liquid to determine the first-order melting transition and its jump in entropy. 
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The pancake-vortex lattice in layered superconductors 
defines a tunable soft matter system with astonishing 
properties Among them, the thermodynamic phase 
transition of vortex-lattice melting and its first-order 
character is now experimentally well established ||J^ , but 
questions remain as to which correlations are lost at the 
transition [Q. Theoretically, the position of the melt- 
ing line can be estimated with a Lindemann criterion 
but a more detailed description of melting is re- 
quired to determine the characteristics of the transition. 
The challenge in defining a theoretical scheme describing 
vortex-lattice melting follows from the complexity of the 
vortex system in real superconductors combined with the 
general lack of exact theories of melting. 

In a moderately anisotropic material, such as 
YBa2Cu307, the vortex crystal melts to a line liquid 
and numerical simulations have treated this in detail . 
In Bi2Sr2CaCu208 (BSCCO), however, the couphng be- 
tween layers is so weak that the layered structure (with 
spacing d) plays a crucial role, and the vortex matter acts 
as a collection of interacting two-dimensional (2D) vor- 
tices, or pancake vortices. Rather than using numerical 
simulations |^,|| , we describe here a novel analytic treat- 
ment to track the melting line through the B-T phase di- 
agram in the extreme anisotropic limit of zero Josephson 
coupling between layers. In this limit the 3D pancake- 
vortex lattice (PVL) remains stable at low temperatures 
due to an attractive electromagnetic interaction between 
pancake vortices in different layers, with range X ^ d 
(A is the in- plane penetration depth). Changing the 
magnetic field tunes the relative importance of this at- 
tractive interlayer interaction and the long-range repul- 
sion between vortices in the same layer. At high fields 
B ^ B\ — <I>o/A^, the in-plane interactions dominate 
and the 3D lattice melts to independent 2D liquids (a 
pancake- vortex gas) close to the 2D melting temperature 
« eod/70 0] [where Eq = {^q/AttX)^]. At lower 
fields the interlayer attraction stabilizes the lattice and 
increases the melting temperature. In the low-field limit 
of weakly-coupled ID stacks, the crystal melts below the 



evaporation transition of an isolated stack of pancake vor- 
tices 1 11 located at the Berezinskii-Kosterlitz-Thouless 
(BKT) vortex unbinding transition of an isolated 
layer at Tbkt = eod/2 - 35T^°. The field regime B - Ba 
where the melting line interpolates between the above 
limits then spans a factor 35 in "reduced" temperature 
T/egd (in real superconductors Sq vanishes as Tf. is ap- 
proached, and the real temperature ratio T^^t/T^ will 
be smaller). Before reaching zero field, the melting line is 
cut by a competing low-field reentrant transition to a di- 
lute liquid of stacks with exponentially-weak interactions 
||l|. A Lindemann analysis ||^ tells us that at T ~ eod/2 
reentrant melting occurs at a field below 10~^B\. 

Ignoring reentrance, we have a 3D melting line that 
interpolates between a ID stack evaporation (exactly de- 
scribed by 2D BKT theory) and a 2D melting transition 
usually described by a BKT-type mechanism of disloca- 
tion pair unbinding. Both limits are well described by a 
self-consistent approximation and we here generalise this 
to all magnetic fields. Our self-consistent method re- 
lies on the long range of the inter-layer attractions; each 
pancake vortex feels the attractive force of pancake vor- 
tices in ~ \/d other layers. Therefore the fluctuations 
in pancake-vortex positions may be averaged, leading to 
an accurate "mean-field" approach where the 2D lattice 
in one layer sits in a substrate potential due to the at- 
traction of the vortex stacks in all other layers. With this 
substrate model we calculate the fluctuations of individual 
pancake vortices, which in turn smears the substrate po- 
tential, and we solve self-consistently. The upper bound 
in temperature to a self-consistent solution leads to an in- 
stability line which we calculate in this paper. We then 
determine the melting line by comparing the free energy 
of the 3D PVL to the free energy of a collection of 2D 
liquids, using numerical results for the 2D system px[ |. 

The evaporation at Tqkt of a single stack of pancake 
vortices occurs because each element is only logarithmi- 
cally bound to the stack: A pancake vortex in a lay- 
ered superconductor generates supercurrents within each 
layer, resulting in a pairwise interaction energy |l3]. 
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where n is the number of layers separating the two vor- 
tices, and R is the in-plane distance. The form of 
this interaction for different limits is well documented 
0: the in-plane repulsion is Vb(i?) = —2eod ln(i?/L) 
(where L is the system-size cut off) and the out-of- 
plane attraction has the large R limit, Vn^o{R) — 
eorf(rf/A)e-"'^/^ln(i?/i). This implies that the energy 
to pull a single pancake vortex (of core-size ^) from a 
straight stack is 2eod ln(i?/^) when R ^ X and the en- 
tropy will unbind the pancake vortices above Tbkt • 
This stack evaporation is easily reproduced within a 
self-consistent substrate model Here, each pancake 
vortex is subject to a quadratic potential, but with a 
strength chosen to match the thermal average ((. . .)) of 
the curvature in the real potential T4(uo) = ^asUo^, 
where u„ is the n-th pancake vortex displacement, and 
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We ignore correlations in the pancake vortex fluctua- 
tions and use the identity for Gaussian fluctuations that 
(exp[-iK • (u„ - uo)]) = exp {-K'^{u'^) /2), to give 
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where Vn{K) = / Re~^^ '^Vn{R). The equipartition 
theorem for a harmonic potential, (w^) — 2T /ag, allows 
us to solve Eq. (^) self-consistently: for large displace- 
ments (u^) ^ }? the limiting form is as = eod/{{u^) + 
2A2), which has the solution (u^) = 2X'^/[l - {2T/eod)], 
diverging at the evaporation temperature Tbkt = Sod/2. 

We now extend this self-consistent analysis to the 
full 3D system at finite fields. We consider the full 
2D fluctuations of the crystal in each layer, sitting on 
a substrate due to the stacked vortex crystals in the 
other layers. Before deriving this in detail, we give 
a quick-and-dirty derivation of evaporation at small 
fields. Close enough to Tbkt the instability occurs 
when <C (u^) ^ Oq, for a vortex density n„ — 
^o/B = 2/^/3aQ. In this limit the substrate poten- 
tial picks up a negative background contribution (see 
below), as ~ SQd[—2Tmy + l/{{u'^)+2X'^)]. Inserting 
this to the equipartition result gives the quadratic equa- 
tion in (^2) = (2r/£od)[2A2 + (u^) + 27Tn,{uY], 
which only has solutions below a temperature given by 
(1 — Tbkt/T')^ — IGTrn^jA^ = 0, and the instability line 
approaches the zero-field transition in the form 
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Note also that this instability occurs when the fluctu- 
ations reach the condition (u^) ~ cqA. This contrasts 



with the often used Lindemann criterion for melting at 
(m^) = c^Og and corresponds to a field dependent Lin- 
demann number cl '■^ (B/Bx)^^'^ [see ^ where a field- 
dependent cl was also found]. 

A precise treatment that can be used at all fields 
must include the elastic distortions of the lattice within 
each layer. Within the self-consistent harmonic approx- 
imation (SCHA) plus substrate model the average en- 
ergy cost for these distortions is given by a quadratic 
form, integrated over all 2D modes in the Brillouin zone, 
Fh[uO] = (I/Stt^) J^^d^Ku"^{K)<!>'^{K)u°{~K) where 
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with and the real and reciprocal lattice vec- 
tors and /y(Q) = QiQj(47r/Q2)e-Q'("')/2. The first 
two terms in (^) are the 2D-elasticity and the last is 
the contribution from the substrate. Again, we have ig- 
nored correlations in displacements (u™(R^)u"(R,y)) = 
<^^i/'^mn<5y (u^)/2 in the last line. Note that the — 
substrate term is cancelled by the = part of the sec- 
ond 2D-elastic term. This is a rcficction of the long range 
divergences in the 2D system with log (Coulomb) inter- 
actions, which do not exist for the 3D system of stacks 
where the circulation currents are screened beyond A. 

The elastic matrix decomposes to transverse 
{5ij - KiKj)^^^K) = cqq{K)K^ + UyUs and longitudi- 
nal KiKj^^^ (K) = cii(K)if2 -f nyUs projections, where 
C66 and cii are the dispersive shear and compression 
moduli, and as is the substrate strength. 
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The (u^) limit of cge and cn recovers the usual 
form for the elasticity of a 2D vortex lattice , with a 
shear modulus Cgg = n„eo(^/4 and a diverging compres- 
sion modulus c'j'i(i^) = 47reorff^^/^^ at small K. A finite 
(m^) softens these moduli, although the diverging com- 
pression modulus remains. To solve self-consistently for 
(m^^ we use the equipartition result with the substrate. 
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where in the last line we have used the small K limits for 
cge and cn and the circular-Brillouin-zone approximation 
with ~ V^/dQ- Note how the substrate potential 
cuts off the log divergence for the 2D shear fluctuations. 
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The self-consistent Eqs. (|^) and determine (m^). 
Above a temperature T„ there are no solutions, giving an 
upper bound to the melting transition. However, in 
it was shown that the SCHA underestimates the degree of 
anharmonic thermal softening of a 2D lattice because of 
the neglect of odd terms in the anharmonicity. Such cu- 
bic anharmonicities are included in the two-vertex-SCHA 
[ p5[ giving results which compare well to numerical simu- 
lations. In Fig. we show the instability line in the B-T 
phase diagram using the two-vertex-SCHA. The line in- 
terpolates between Tbkt at small fields and at high 
fields. Also shown is the significant thermal softening of 
both cgg and as at S = B\ with the resulting anharmonic 
contributions (beyond linear in T) for (u^). 
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FIG. 1. The instability hne for the PVL in the B-T 
plane calculated with the two-vertex-SCHA. The line goes 
asymptotically to the 2D melting temperature at high fields, 
and ends at Tbkt at zero field. The left inset shows the 
two-vertex-SCHA results for the shear modulus cee, the sub- 
strate strength Os, and the pancake fluctuation width (u^) 
with increasing temperature at _B = Bx- The right inset 
is a low-field zoom of the instability line showing the result 
(dashed) of the simple SCHA scheme for comparison: the two 
lines coincide at very low fields when the substrate dominates 
over the 2D elasticity, but the SCHA overestimates the stabil- 
ity of the lattice at higher fields. We do not include here the 
intrinsic temperature dependence in real superconductors of 
the penetration depth and the energy scale eo (but see Fig. 2) . 

The above instability line only marks the upper tem- 
perature limit of the solid phase. A true first-order tran- 
sition occurs when the free energies of two phases cross. 
For melting transitions it is often difficult to accurately 
determine the free energy of the liquid. We can make 
progress for PVL melting, as the liquid state behaves 
to a good approximation as uncoupled 2D liquids in each 
layer. The free energy of a 2D liquid with log-interactions 
can be extracted from numerical simulations ]lO| . A 
crucial ingredient is that the free energy is known ex- 
actly at one special temperature T = e^d. For the 
solid we calculate the free energy of a 2D crystal on a 



self-consistent substrate. We must be careful to take 
the same normalization in both phases: as in jl^ we 
normalize with respect to the ideal gas of N pancakes, 
Z{H = 0) = 1, and define the partition function. 
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fixing the free energy, measured from the ideal gas, 
as F = -TlnZ. We write H ^ H^u + Hs where 
= {N/2)[j:^^^oVo{R^)-ny J d^RVoiR)] and i/^ = 
iN/2) [Er,,„^o VniRf^) + nvj d^RVoiR)] are the 2D in- 
teraction term and substrate energy respectively . 

In the solid, the free energy of harmonic fluctuations 
is straightforward to calculate. The right-hand side of 
the inequality F < -\- {H — H^)\i is minimized by 
the SCHA (see where is defined by (|). The 
harmonic free energy of 2D fluctuations is 
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The correction part of the variational free energy has con- 
tributions from and Hg. Ignoring the small anhar- 
monic part of the 2D energy, the difference {H^d — i?^D)h 
is just the ground state energy (which is not included in 
H^j^). In the energy for a 2D lattice of log- interacting 
particles is found to be E^^ = 0.7A9Neod. The substrate 
energy correction. Eg = {Hs — Hl^)h, is 
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and the sum F^ + E^^ + Eg is the variational free energy. 

In the liquid phase, {Hg) is very small and we should 
find the free energy of a 2D liquid. The internal energy 
U{V) (with F = 2eod/T) is found from simulations Q, 



C/(r) = (-0.751 + 0.880 F 



-0.74 



0.209 F 



-l.7^ 



Neod. (11) 



From the relation F + TdrF, the excess free energy 
can be written as FFiiq(r) = FoFiiq(ro) + Jp^ J7(F')rfF'. 
We know the value at Fo = 2, where Fiiq(2) = °0.081eodN 
|0, so we can integrate (|0) to give (with t ~ T/e^d), 
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0.092ii-^)iVeod. (12) 



Using these results we can plot the free energy of both 
phases and see where they cross; a typical example at 
B — Bx is shown in the inset of Fig. ^. Calculating the 
crossing point numerically at different fields, we find a 
melting line (see Fig. ||) just below the stability limit of 
the lattice that approaches zero field at T w 0.47eoc'- 

The jump in slope S — ^OtF gives a latent heat 
TAS" = AC/. In Fig. |^ we plot the entropy jump per pan- 
cake vortex as a function of the transition temperature. 
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At high fields (small T) it approaches the value As„ « 
0.4fcB, consistent with 2D simulations and simple es- 
timates At low fields (T ~ O.beod) the latent heat 
appears to weakly diverge as T — > Tbkt- We understand 
this as the energy of the liquid is roughly constant, U{T m 
Tbkt) ~ —Q-2l%Neod, while the energy of the soHd is 
dominated by the substrate term, Eg ~ iVeoc?ln(A/ao) so 
that the latent heat is AJ7 ~ —Es- This gives an entropy 
jump per vortex pancake Asu = AU/NT ^ ln{B\/B). 
This is of the same form as the entropy difference be- 
tween an ideal gas, Asf,'^^ « A:Bln(aQ/^^) and the (re- 
duced phase-space) solid As^°' w fee In ((u^)/^^) when 
(m^) ^ qqX (as found above for the low-field instability). 
We do not include here the T dependence of A in real 
systems that gives extra terms in the latent heat [^ . 
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FIG. 2. The full line in the upper graph shows the melt- 
ing line Bm{T) of the pancake vortex lattice as calculated by 
comparing the free energies of the solid and liquid phases. 
The dashed line is the stability limit Bu{T) of the PVL as 
shown in Fig. ^. The inset shows the free energy comparison 
at B = B\. The lower graph gives the entropy jump per pan- 
cake vortex Asv Also shown are the real scales of T and B 
assuming A(r) = A(0)/[1 - {T/Tcof]2, with A(0) ^ 2000A, 
d « I5A and Tco ~ 100 K. The low-field reentrant melting 
line, not shown here, will cut Bm{T) below 10~^i3A. 

Previously, melting of the magnetically coupled PVL 
has been analyzed numerically and via density func- 
tional theory (DFT) of the liquid phase [|o|. The early 
simulations in find evidence of melting at T « O.OOeod 
when B/Bx ~ 0.2 (close to our melting line). The DFT 
gives the stability limit of the liquid and provides results 
consistent with ours at fields above 0.5i?A- At lower fields 
the DFT gives a liquid instability above our melting line, 
a discrepancy which requires further study. 

To compare to real superconductors our units of field 
B\ and of temperature eod must be scaled due to the 
intrinsic variation in A(T) (diverging at a temperature 



T°). Doing this, we find that our melting curve of Fig. ^ 
lies below experimental melting lines |2pll| for reason- 
able choices of A(T); this is because the Josephson cou- 
pling, neglected in this paper, becomes important as T° 
is approached Q and stiffens the vortex lattice. Our re- 
sults may be recovered in experiments by suppressing the 
Josephson coupling with a strong in-plane magnetic field. 
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